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PREFACE 
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Mechanics,"  Task  2307P102,  "Failure  Prediction  in  Metals."  Professor  J. 

N.  Yang  was  engaged  as  a  visiting  scientist  from  the  School  of  Engineering 
and  Applied  Science,  The  George  Washington  University,  Washington,  D.C. 
20052. 

This  report  covers  research  conducted  during  the  period  February 
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SECTION  I 
INTRODUCTION 

Service  loading  spectra  experienced  by  gas  turbine  engine  disks 
include  not  only  fatigue  cyclic  loadings  but  also  sustained  loads  for 
various  periods  of  time.  Crack  growth  damage  accumulation  for  engine 
materials  occurs  not  only  under  load  reversals  but  also  under  sustained 
loads  combined  with  high  temperatures  and  corrosive  environments.  While 
the  crack  growth  behavior  of  engine  materials  due  to  the  interactions 
between  fatigue  and  sustained  loads  has  not  been  understood  well, 
analytical  and  experimental  studies  have  been  performed  under  sustained 
loads  (e.g..  References  1  and  2)  or  fatigue  load  alone  (e.g..  Reference 
3)  at  high  temperatures  based  on  the  fracture  mechanics  approach. 

An  experimental  test  program  for  the  crack  growth  of  IN100,  a 
superalloy  used  in  F100  engine  disks,  under  sustained  loads  at  elevated 
temperature  has  been  conducted  (References  1  and  2).  Test  results  of 
the  crack  growth  rate  using  23  compact  tension  specimens  obtained  in 
References  1  and  2  indicate  a  considerable  statistical  scatter.  As  a 
result,  the  statistical  treatment  for  the  crack  growth  variability  is 
necessary  to  permit  a  rational  life  prediction  of  gas  turbine  engine 
disks,  in  particular  the  effect  of  hold  time  in  the  load  cycle.  A 
fracture  mechanics-based  statistical  model  for  the  crack  growth  damage 
accumulation  of  IN100  under  sustained  load  is  proposed  and  investigated. 

Specimen  crack  length  records  versus  time  were  generated  from  load 
line  displacement  versus  time  data.  These  data  were  smoothed  against 
time  along  with  its  first  derivative  using  a  seven-point  averaging  method 
making  use  of  a  least  squares  fit  to  a  second-degree  polynomial.  By 
using  the  above,  together  with  experimentally  determined  compliance 
versus  time  for  each  specimen,  an  "effective"  modulus  of  elasticity  for 
each  specimen  was  determined  and  a  smoothed  effective  crack  length  versus 
time  record  was  obtained.  This  accounted,  in  part,  for  crack  front  tunnel¬ 
ing  which  in  IN100,  when  under  sustained  load  at  elevated  temperature, 
appears  to  be  significant.  The  details  of  data  reduction  are  discussed 
in  Reference  2. 
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The  crack  growth  rate  data  alone  are  utilized  to  calibrate  the 
fracture  mechanics  parameters  and  the  statistical  properties  of  the  model 
using  the  method  of  maximum  likelihood.  A  correlation  study  is  carried 
out  for  the  crack  length  as  a  function  of  time.  From  the  statistical 
model,  the  distributions  of  the  crack  growth  rate,  the  crack  growth  life 
to  reach  any  given  crack  size,  and  the  crack  size  at  any  time,  have  been 
derived.  It  is  demonstrated  that  the  correlation  between  the  experimental 
test  results  and  the  statistical  model  is  very  reasonable. 
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SECTION  II 
TECHNICAL  APPROACH 

Several  mathematical  models  were  investigated  for  the  description 
of  the  crack  growth  behavior  of  IN100  under  sustained  loads  at  elevated 
temperatures  based  on  the  principle  of  fracture  mechanics  (References  1 
and  2).  These  models  can  be  expressed  in  a  general  form 

da(t)/dt  =  l(K,G,J,C*)  (1) 

in  which  a(t)=crack  length  at  time  t,  L=a  general  non-negative  function, 
K=stress  intensity  factor,  G=strain  energy  release  rate,  J=J-integral , 
and  C*=a  line  integral  related  to  the  rate  of  change  of  potential  energy 
per  unit  of  crack  growth.  A  special  form  of  L  in  Equation  1  was  observed 
to  be  reasonable  for  the  crack  growth  rate  data  of  IN100  CT  specimens. 

Mil  =  QKb  (2) 

However,  even  under  well-controlled  laboratory  conditions,  crack 
growth  rate  data  exhibit  a  considerable  statistical  variability.  In  the 
life  prediction  of  engine  components  such  a  statistical  scatter  should 
be  taken  into  account  rationally;  thus  a  statistical  model  based  on 
Equation  2  is  needed.  We  propose  to  randomize  Equation  2  as  follows; 

M^-  =  QKbX(K)  (3) 

in  which  X(K)  is  a  non-negative  random  function  (or  process)  of  K  taking 
values  around  1.0  to  account  for  the  random  variation  of  the  crack  growth 
rate.  Hence,  the  deterministic  fracture  mechanics  model  of  Equation  2 
represents  the  mean  crack  growth  behavior  whereas  the  statistical 
properties  of  the  crack  growth  rate  are  taken  care  of  by  the  random 
function  X (K ) - 

Two  extreme  cases  of  the  random  function  X ( K )  should  be  mentioned. 

At  one  extreme  X(K)  is  totally  uncorrelated  at  any  two  different  values 
of  K.  Based  on  the  central  limit  theorem,  the  statistical  dispersion  of 
the  crack  length  a(t),  after  integrating  Equation  3,  is  the  smallest 
among  the  class  of  random  functions  X(K)  (References  4,  5).  Hence,  it 
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is  unconservative  for  engineering  analysis  and  design  purposes.  At 
another  extreme,  the  random  function  X(K)  is  completely  correlated  at 
all  values  of  K,  indicating  that  X(K)  is  a  random  variable  X.  When  X(K) 
is  a  random  variable,  a  specimen  with  a  higher  (or  lower)  crack  growth 
rate  will  maintain  it  throughout  the  entire  life  time.  The  totally  cor¬ 
related  random  function  X(K)=X  results  in  the  largest  statistical  dis¬ 
persion  of  the  crack  length  a(t).  Hence,  it  is  conservative  for  life 
prediction  and  engineering  application.  It  should  be  mentioned  that 
instead  of  a  random  function  of  K,  a  random  function  of  time  t,  i.e., 

X(t) ,  can  also  be  used  (Reference  6). 

When  X(K)  is  a  random  function,  the  crack  length  a(t)  is  also  a 
random  function  of  time  t  as  a  consequence  of  Equation  3.  The  solution 
for  the  statistical  properties  of  a(t)  is  rather  involved  and  it  may  be 
beyond  the  means  of  current  engineering  discipline  (Reference  6). 

Because  of  its  simplicity  and  conservative  nature  for  crack  propagation 
prediction,  we  shall  explore  the  extreme  case  in  which  X(K)=X  is  a  random 
variable,  i.e., 

MU  =  QKbx  (4) 

in  which  X  is  a  positive  random  variable  taking  values  around  unity. 

Taking  the  logarithm  of  both  sides  of  Equation  4,  one  obtains 


Y  =  bU  +  q  +  Z  (5) 

where 


Y  =  log  Mil  ,  U  =  log  K  1 
dt  I 

(6) 

Z  =  log  X  ,  q  =  Log  Q  J 

where  Z  is  a  random  variable  taking  values  around  zero. 


Now  Z  is  assumed  to  be  a  normal  random  variable  with  zero  mean  and 
the  standard  deviation  a z.  Then  it  follows  from  Equation  5  that  the  log 
crack  growth  rate  Y  =  log  da(t)/dt  is  also  a  normal  random  variable  with 
the  mean  value  uy  and  standard  deviation  a y  given  by 
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uy  =  bU  +  q  (7) 

lly  =  '’z  <8> 

The  parameters  b  and  q  (or  Q)  as  well  as  the  standard  deviation 
(or  o  )  can  be  estimated  from  the  test  results  of  the  crack  growth  rate 
versus  stress  intensity  factor  using  Equation  5  and  the  method  of 
maximum  likelihood.  Since  Y  and  Z  are  normal  random  variables  and 
Equation  5  is  linear,  the  method  of  maximum  likelihood  is  identical  to 
the  linear  regression  analysis  and  the  method  of  least  squares. 

Test  results  of  the  crack  growth  rate  versus  stress  intensity 
factor  for  I N 100  CT  specimens  at  1350°F  have  been  generated  for  five 
different  thicknesses.  The  results  are  presented  in  Figures  1(a) 
through  1(e).  The  method  of  maximum  likelihood  has  been  used  to 
estimate  t'he  values  of  b,  Q,  and  oy=az  for  five  test  groups  and  the 
results  are  shown  in  Table  1  on  page  9.  The  straight  lines  in  Figures 

1(a)  through  1(e)  represent  the  estimated  mean  values  u y =  bll  +  q 
(Equation  7). 

Each  datum  in  Figures  1(a)  through  1(e)  is  calculated  from  a  one- 
to-one  corresponding  crack  growth-versus-time  and  specimen  load  line 
displacement-versus-time  measurement  using  a  seven-point  averaging 
method  to  smooth  the  time  records.  Crack  lengths  were  computed  from 
experimental  compliance  based  on  an  exact  solution  for  load  line 
displacements  (Reference  2). 

The  initial  crack  growth  rates  for  the  23  specimens  were  observed 
to  start  from  an  initially  high  rate  and  then  decrease  to  a  level  which 
corresponded  to  the  onset  of  stable  creep  crack  growth.  This  decrease 
in  the  initial  crack  growth  rates  was  due  to  an  initially  decreasing 
load  line  displacement  rate,  upon  which  the  calculations  of  crack 
lengths  rested.  These  initially  decreasing  rates  were  excluded  from  the 
analysis  performed  herein. 

To  show  the  validity  of  the  assumption  that  Z  follows  the  normal 
distribution,  sample  values  of  Z  denoted  by  Zy  are  computed  from  the 
test  results  of  the  log  crack  growth  rate  Y  -  log  da(t)/dt  versus  the 
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Log  K.  Ksi.in**l  2 

ure  1(a).  Log  Crack  Growth  Rate  Versus  Log  Stress  Intensity  Factor 
for  Various  Thickness  Specimen  Group;  7/32  Inches. 


0 . 5 


l 


1 1 /32  Inches 


Test  Condition  No.  2 


Figure  1(b), 


Log  K.  Ksi.in*M  2 

Log  Crack  Growth  Rate  Versus  Log  Stress  Intensity  Factor 
for  Various  Thickness  Specimen  Group;  11/32  Inches. 
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1.3  1.4  1.3  1.6  1.7  1.8  1.9  2.8 

Log  K.  Ksi.in**l  2 

Figure  1(c).  Log  Crack  Growth  Rate  Versus  Log  Stress  Intensity  Factor 
for  Various  Thickness  Specimen  Group;  15/32  Inches. 


Log  K,  Ksi.in**  1  2 

Figure  1(d).  Log  Crack  Growth  Rate  Versus  Log  Stress  Intensity  Factor 
for  Various  Thickness  Specimen  Group;  19/32  Inches. 
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log  stress  intensity  factor  U  =  log  K,  denoted  by  (y.,  u-)>  using 

J  J 

Equation  5 


2j  =  yj  '  buj  "  q  for  j=1»2*--**n  (9) 

in  which  b  and  q  have  been  estimated  by  the  method  of  maximum  likelihood 
and  n  is  the  total  number  of  test  data. 

Sample  data  of  Zj  ( j=l,2, . . . ,n)  obtained  from  Figure  1(b)  are 
plotted  on  normal  probability  paper  in  Figure  2  as  circles  along  with  a 
straight  line  representing  the  estimated  normal  distribution  of  Z  (with 
zero  mean  and  standard  deviation  determined  previously,  see  Table  1). 

A  linear  scale  is  used  in  Figure  2  in  which  the  sample  data  z-  are 

J 

arranged  in  ascending  order,  i.e.,  z1<z,<z,. . ,<z  ,  and  the  ordinate 

i  i  c.  o  n  ^ 

corresponding  to  z-  is  given  by  j/(n+l) ]  with  $'i(  )  being  the 

J 

inversed  standardized  normal  distribution  function.  It  is  observed  from 
Figure  2  that  the  sample  values  of  Z  are  scattered  around  the  straight 
line  without  a  nonlinear  trend,  indicating  that  the  normal  distribution 
is  acceptable. 


TABLE  1 


MAXIMUM  LIKELIHOOD  ESTIMATE  OF  b,  Q,  STANDARD  DEVIATION  o  AND 
COEFFICIENT  OF  VARIATION  OF  da/dt  y 


Thickness 

(Inch) 

b 

Q 

a  =o 

y  z 

Coef.  of 
Variation 

No.  of 
Points 

No.  of 
Specimens 

7/32 

3.1015 

4.8563xl0'7 

0.13325 

31.42% 

233 

4 

11/32 

2.1626 

2.7349xl0'5 

0.06796 

15.75% 

227 

5 

15/32 

2.4995 

8.6942xl0"6 

0.09550 

22.30% 

156 

4 

19/32 

2.5255 

8.0931xl0'6 

0.12870 

30.30% 

185 

5 

23/32 

2.4220 

1.2594xl0"5 

0.05740 

13.28% 

132 

4 

The  Kolmogorov-Smirnov  test  for  goodness  of  fit  has  been  performed 

to  determine  the  observed  K-S  statistic  D  .  It  is  shown  that  the  normal 

n 

distribution  is  acceptable  at  least  at  a  10%  level  of  significance. 

Since  Z=log  X  and  Y=log  da(t)/dt  are  normal  random  variables,  the 
crack  growth  rate  G=da(t)/dt  follows  the  lognormal  distribution.  The 
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Figure  2.  Normal  Probability  Plot  for  Z  for  11/32  Inches  Thick  Specimen 
Group. 
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coefficient  of  variation,  V,  of  da(t)/dt  is  related  to  the  standard 
deviation  oy =a2  through 

V  =[e(°y1nl0)2  -  1]  1/2  (10) 

Test  results  of  the  coefficient  of  variation  V  for  the  crack  growth  rate 
G=da(t)/dt  for  five  test  groups  are  presented  in  Table  1. 

It  is  interesting  to  mention  that  the  coefficient  of  variation,  V, 
which  is  a  measure  of  the  statistical  dispersion  of  the  crack  growth 
rate  for  IN100  CT  specimens  under  fatigue  cyclic  loadings  at  elevated 
temperatures  reported  in  Reference  4,  is  larger  than  that  of  the  same 
specimen  under  sustained  loads  obtained  in  this  report. 

The  distribution  of  Z  and  Y=log[da( t)/dt]  are  both  normal  with  the 
same  standard  deviation  az=oy  that  has  been  determined  previously.  Let 
zy  be  the  y  precentile  of  Z,  i.e., 

y%  =  P[Z>7y]  =  l  -  o(zy/o2)  (11) 

or  inversely, 

*Y  =  oz  ^(l  -  Y%)  (12) 

The  y  percentile  of  the  log  crack  growth  rate  Y,  denoted  by  yy  follows 
from  Equation  5  as 

y  =  bU  +  q  +  z  (13) 

Y  Y 

in  which  zy  is  given  by  Equation  12. 

Therefore,  by  varying  the  value  of  Y,  one  obtains  the  distribution 
of  the  log  crack  growth  rate  in  terms  of  percentiles.  The  results  are 
shown  in  Figure  3  for  the  second  test  group.  As  an  example,  the  crack 
growth  rate  path  associated  with  y=10  indicates  that  10%  of  specimens 
will  have  a  growth  rate  faster  than  that  shown  by  the  curve. 

The  y  percentile  of  the  random  variable  X,  denoted  by  XY,  is 
computed  from  Equation  6  as 

XY  =  (10)Zy 

in  which  z  is  given  by  Equation  12. 


(15) 
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Figure  3.  Percentiles  of  Log  Crack  Growth  Rate  Y  as  a  Function  of  Log 
Stress  Intensity  Factor  X  for  11/32  Inches  Thick  Specimen 
Group  (Theoretical  Model J".' 
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The  y  percentile  of  the  crack  length  a(t)  at  time  t,  denoted  by 
ay(t),  is  obtained  by  substituting  Equation  15  into  Equation  4 


da  (t) 
dt 


QKbX 

y 


(16) 


in  which  the  stress  intensity  factor  K  of  the  compact  tension  specimen 
is  expressed  as 

(17) 


where  B=thickness  of  the  specimen,  W=width  of  the  specimen,  P=applied 
load  a  id 


f  [ay(t)] 


(0.866  +  4.64a  -  13.32a2 


(1-a) 


3/2 


14.72a3  -  5.6a4)  (18) 


a  (t) 
a  =  _ 

W 


(19) 


Thus,  Equation  16  can  be  integrated  numerically  over  particular 
limits  to  obtain  a  set  of  crack  length  ay(t)  versus  time  t  for  different 
Y  percentiles  as  shown  in  Figure  4.  The  integrations  were  performed  for 
an  ASTM  compact  tension  specimen  with  thickness  B  =  11/32  inch  (8.73  mm), 
width  U  =  2.5  inches  (63.5  mm),  maximum  load  P  =  3.17  kips  (14.1 kN ) , 
initial  crack  length  a( 0)  =  0.68  inch  (17.3  mm)  and  final  crack  length 
a(t)  =  1.25  inch  (31.75  mm).  It  should  be  mentioned  that  the  numerical 
integration  of  Equation  16  for  each  y  percentile  of  the  crack  length, 
ay(t),  is  deterministic  and  straightf orward. 

In  the  life  prediction  of  engine  components,  two  statistical  dis¬ 
tributions  are  of  practical  importance;  (1)  the  distribution  function, 

Ff( a) ( t the  time,  T(a),  to  reach  any  given  crack  length,  a,  and  (2) 
the  distribution  function,  Fa^(u),  the  crack  length,  a(t),  at  any 
time  t.  It  has  been  shown  in  Reference  4  that  both  distributions  are 
related  and  one  can  be  derived  from  the  other. 


Crack  Length,  in. 
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1 . 05 
0 . 95 
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0 . 75 
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Time,  hrs. 

Figure  4.  Percentiles  of  Crack  Length  a{t}  as  Function  of  Time  t 

(Theoretical  Model)  for  11/32  Inches  Thick  Specimen  Group. 
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Since  Figure  4  represents  the  distribution  of  the  crack  length  as  a 
function  of  time  t,  it  contains  all  the  information  needed  to  determine 
the  distributions  mentioned  previously.  For  instance,  by  drawing  a 
horizontal  line  in  Figure  4  through  a  crack  length  of  interest,  the  dis¬ 
tribution  of  time  to  reach  that  crack  length  is  obtained.  Likewise, 
drawing  a  vertical  line  in  Figure  4  through  a  given  time  t,  one  obtains 
the  distribution  of  the  crack  length  at  that  time  instant.  The  complement 
F*a(t)(u)  of  the  distribution  function  Fa^(u)  of  the  crack  length, 
i.e.,  fr*a(t)Cu)  =  l=fra(t)^u^*  is  the  Probability  that  the  crack  length  at 
time  t  will  exceed  certain  length  u.  Hence,  F*a(t)^u)  referred  to  as 
the  crack  exceedance  curve. 

The  distribution  function  FT^(t)  tor  the  time  to  reach  the  crack 
lengths  a=0.75,  0.85,  0.95,  and  1.25  inches  (19.05,  21.6,  24.1  and  31.75 
mm,  respectively)  are  obtained  from  Figure  4  and  displayed  in  Figure  5 
as  solid  curves.  The  crack  exceedance  curve  F*a(t)^  at  hours  is 

shown  in  Figure  6  as  a  solid  curve. 

It  should  be  emphasized  that  the  numerical  integration  of  Equation 
16  for  each  r  percentile  of  the  crack  length,  ay(t),  is  deterministic 
and  straightforward.  As  a  result,  it  is  very  easy  to  construct  the  dis¬ 
tribution  of  the  crack  length  as  a  function  of  time  given  by  Figure  4. 
Thus,  the  present  statistical  model  is  very  simple  for  practical 
appl ications. 
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Figure  5.  Distribution  of  Time  to  Reach  Given  Crack  Lengths  a  =  0.75, 
0.85,  0.95,  and  1.25  Inches  for  11/32  Inches  Thick  Specimen 
Group  (Solid  Curves  are  Theoretical  Model). 


Figure  6.  Probability  of  Crack  Exceedance  at  2.5  Hours  for  11/32  Inches 
Thick  Specimen  Group  (Solid  Curve  is  Theoretical  Model). 
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SECTION  III 

CORRELATION  WITH  TEST  RESULTS 

A  problem  in  the  correlation  study  is  that  the  number  of  specimens 
in  each  thickness  group  is  very  small  (Table  I)  and  even  with  such  a 
small  sample  size  each  specimen  is  usually  tested  using  different  initial 
flaw  length,  a(0),  final  flaw  length,  ap,  and  applied  load  P.  The  load¬ 
ing  conditions  for  each  specimen  were  given  in  Reference  2.  Thus  the 
test  results  of  the  crack  length  a(t)  are  not  homogeneous. 

Fortunately,  the  present  statistical  model  utilizes  only  the  crack 
growth  rate  versus  the  stress  intensity  data  (Figure  1)  to  calibrate  all 
the  fracture  mechanics  parameters  without  using  the  test  results  of  the 
crack  length  a ( t )  versus  time  t.  To  study  the  correlation  between  the 
statistical  model,  and  the  test  results,  in  terms  of  the  crack  growth 
damage  accumulation,  a(t),  versus  time,  t,  including  the  distribution  of 
time  to  reach  any  given  crack  length,  and  the  distrioution  of  the  crack 
length  at  any  time  instant  t,  homogeneous  test  conditions  were  assumed 
for  each  thickness  group  as  shown  in  Table  2.  The  applied  loads  P  used 
in  Table  2  were  chosen  to  avoid  excessive  extrapolation  into  the  region 
of  the  stress  intensity  factor  K  where  no  actual  data  existed  (Figures 
l(a)-l(e)). 

TABLE  2 

ASSUMED  HOMOGENEOUS  TEST  ENVIRONMENTS  FOR 
STANDARD  COMPACT  TENSION  TEST  SPECIMENS 


Thickness 

a  ( 0 ) 

dF 

P 

(in.) 

(in. ) 

(in.) 

(kips) 

7/32 

0.66 

1.25 

2.28 

11/32 

0.68 

1.25 

3.17 

15/32 

0.68 

1.25 

4.39 

19/32 

0.68 

1.25 

4.08 

23/32 

0.68 

1.25 

5.89 

Note:  a( 0)  =  Initial  Flaw  Length,  ap  =  Final  Flaw  Length, 
P  =  Appl ied  Load 
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Test  results  of  the  crack  growth  rate  for  each  specimen  are  best 
fitted  by  Equation  2  using  the  method  of  least  squares  tc  estimate  the 
values  of  o  and  Q.  Sample  values  of  b  and  Q  for  each  specimen  are 
presented  in  Table  3  for  various  thickness  group.  The  best  fitted  crack 
growth  rate  's  shown  in  Figure  7  as  a  straight  line  for  the  first  speci¬ 
men  in  eu^n  thickness  group,  along  with  the  test  data.  It  should  be 
mentioned  that  the  crack  growth  curves,  a(t)  versus  time  t,  obtained  by 
integrating  the  test  data  of  Figure  7,  are  indistinguishable  from  those 
determined  by  integrating  the  best  fitted  straight  line.  In  other  words, 
the  resulting  crack  length  a(t)  versus  time  t  obtained  by  integrating 
Equation  2  and  using  the  values  of  b  and  Q  given  in  Table  3  are  essentially 
identical  to  that  obtained  by  integrating  the  actual  crack  growth  rate 
data.  This  justifies  the  correlation  study  of  the  crack  length  using 
the  assumed  homogeneous  test  conditions  given  in  Table  2. 

TABLE  3 


BEST  FITTEO  b  AND  Q  VALUES  FOR  EACH  SPECIMEN  IN  ALL  THICKNESS  GROUPS 


:imt-n 

No.  b 

Q 

Specimen  No 

b 

Q 

7/32 

Inches  Thickness  Group 

15/32 

Inches  Thickness  Group 

1 

2.6559 

3.4359xl0-6 

1 

2.9257 

1 . 8229x10*6 

2 

3.0959 

3.4359xl0"7 

2 

2.5697 

7 ,6367xl0*6 

3 

2.8626 

1.4920xl0*6 

3 

1.8834 

7.2033x10*6 

4 

2.4962 

5 .7090x10*6 

4 

2.2316 

2. 4520xl0*5 

11/32 

Inches  Thickness  Group 

19/32 

Inches  Thickness  Group 

1 

2.1847 

2.1388xl0*5 

1 

3.7926 

6.3718xl0”8 

2 

1.4166 

5.5907xl0'4 

2 

2.4683 

1.0241x10*6 

3 

2.1912 

2.6118xl0*5 

3 

2.1385 

5.1545x10*6 

4 

2.7300 

3.7005xl0~6 

4 

2.1953 

2.2517x10*6 

5 

2.6969 

2.9i39xl0~6 

5 

2.9969 

1. 4040xl0*6 

23/32 

Inches  Thickness  Group 

1 

1.9285 

8.3606xl0~5 

2 

2.7192 

4.1785xl0~6 

3 

2.6142 

5.9990xl0~6 

4 

2.6729 

4.8845xl0~6 

Log(da/dt).  ln/^r  £  Log(da/dt),  m,  hr 
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jre  7(a).  Best  Fitted  Crack  Growth  Rate  Data  for  the  First  Specimen 
in  Each  Thickness  Group;  7/32  Inches. 


Log  K.  Ksiin**l  2 

Figure  7(b).  Best  Fitted  Crack  Growth  Rate  Data  for  the  First  Specimen 
in  Each  Thickness  Group;  11/32  Inches. 
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Log  K.  Ksi.in**l  2 

jre  7(c).  Best  Fitted  Crack  GroWth  Rate  Data  for  the  First  Specimen 
in  Each  Thickness  Group;  15/32  Inches. 
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Log  K,  Ksi.in**l/2 

Figure  7(d).  Best  Fitted  Crack  Growth  Rate  Data  for  the  First  Specimen 
in  Each  Thickness  Group;  19/32  Inches. 
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Figure  7(e).  Best  Fitted  Crack  Grc 
in  Each  Thickness 
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On  the  basis  of  the  statistical  model,  the  maximum  likelihood  esti¬ 
mate  of  b,  Q  and  obtained  in  Table  1  are  used  to  construct  various 
percentiles  of  the  log  crack  growth  rate  Y= 1 og  da/dt  versus  log  stress 
intensity  factor  U^log  K.  The  results  for  five  thickness  groups  are 
shown  in  Figures  8(a)-12(a)  for  y  =  5,  10,  25,  50,  75,  90  and  95  percen¬ 
tiles.  The  least  square  estimates  of  b  and  Q  for  each  individual  specimen 
given  in  Table  3  are  used  in  conjunction  with  Equation  2  to  construct  Y 
versus  U,  and  the  results  are  displayed  in  Figures  8(b)-12(b).  Figures 
8(b)-12(b)  are  referred  to  as  the  extrapolated  test  results,  since  the 
actual  test  data  do  not  exist  in  the  entire  range  of  K. 

The  crack  growth  rates  shown  in  Figures  8-12  have  been  integrated 
over  the  limits  given  in  Table  2,  the  assumed  homogeneous  test  conditions, 
to  yield  crack  length,  a(t),  versus  time,  t.  The  corresponding  results 
are  depicted  in  Figures  13-17  for  five  thickness  groups  (again  curves  in 
Figures  13( a) -17( a)  represent  <  =  5,  10,  25,  50,  75,  90  and  95  percentiles). 
A  comparison  between  Figures  13(a)~17(a)  and  Figures  13(b)-17(b)  indicates 
that  the  correlation  between  the  statistical  model  and  the  extrapolated 
test  results  is  very  reasonable. 

Based  on  the  statistical  model,  the  distribution  function,  F j ^ a) ( t) , 
for  the  time  T ( a )  to  reach  crack  lengths  0.75,  0.85,  0.95,  and  1.25  inches 
(19.05,  21.6,  24.1  and  31.75  mm,  respectively)  are  obtained  by  drawing 
horizontal  lines  through  appropriate  crack  lengths  in  Figures  13(a)- 
17(a)  and  the  results  are  plotted  in  Figures  18-22  as  solid  curves 
designated,  respecti vely,  by  Curves  1,  2,  3  and  4.  The  corresponding 
distribution  functions  for  the  extrapolated  test  results  are  obtained 
from  Figures  13(b)-17(b)  and  displayed  in  Figures  18-22  as  circles, 
squares,  solid  circles  and  triangles,  respectively.  In  constructing  the 
distribution  function  for  the  extrapolated  test  results  shown  in  Figure 
18  for  instance,  a  horizontal  line  is  drawn  through  the  given  crack  length 
in  Figure  13(b)  to  obtain  four  data  points.  These  four  data  points  are 
ranked  in  an  ascending  order  and  the  ordinate  of  the  ith  data  point  is 
i/(m+l)  where  m=4  is  the  total  number  of  points.  Figures  18-22  demon¬ 
strate  a  good  correlation  between  the  statistical  model  (solid  curves 
indicated  by  1,  2,  3  and  4)  and  the  extrapolated  test  results  (circles, 
squares,  solid  circles  and  traingles). 
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Log  K,  Ksi.in**l/2 

jre  8(a).  Log  Crack  Growth  Rate  Y  Versus  Log  Stress  Intensity  Factor 
for  7/32  Inches  Thickness  Group;  Statistical  Model. 


Log  K,  Ksi.ir»**1/2 

Figure  8(b).  Log  Crack  Growth  Rate  Y  Versus  Log  Stress  Intensity  Factor 
for  7/32  Inches  Thickness  Group;  Extrapolated  Test 
Results. 
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Log  K.  Ksi.in**l/2 

Figure  9(a).  Log  Crack  Growth  Rate  *Y  Versus  Log  Stress  Intensity  Factor 
for  11/32  Inches  Thickness  Group;  Statistical  Model. 


Log  K.  Ksi.in**l/2 

Figure  9(b).  Log  Crack  Growth  Rate  Y  Versus  Log  Stress  Intensity  Factor 
for  11/32  Inches  Thickness  Group;  Extrapolated  Test 
Results. 
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Log  K,  Ksi.in**l  2 

jure  10(a).  Log  Crack  Growth  R^te  Y  Versus  Log  Stress  Intensity 
Factor  for  15/32  Inches  Thickness  Group;  Statistical 
Model . 
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Log  K,  Ksi.in**l/2 

Figure  10(b).  Log  Crack  Growth  Rate  Y  Versus  Log  Stress  Intensity 
Factor  for  15/32  Inches  Thickness  Group;  Extrapolated 
Test  Results. 
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Log  K.  Ksi.in**l  2 

Figure  11(b).  Log  Crack  Growth  Rate  Y  Versus  Log  Stress  Intensity 
Factor  for  19/32  Inches  Thickness  Group;  Extrapolated 
Test  Results. 
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Figure  12(a).  Log  Crack  Growth  Rate  Y  Versus  Log  Stress  Intensity 
Factor  for  23/32  Inches  Thickness  Group;  Statistical 
Model . 


Log  K.  Ksi.in**l  2 

Figure  12(b).  Log  Crack  Growth  Rate  Y  Versus  Log  Stress  Intensity 
Factor  for  23/32  Inches  Thickness  Group;  Extrapolated 
Test  Results. 
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Time,  hrs. 

ire  13(a).  Crack  Length  a(t)  Versus  Time  t  for  7/32  Inches  Thickness 
Group;  Statistical  Model. 


Time.  hrs. 

Figure  13(b).  Crack  Length  a(t)  Versus  Time  t  for  7/32  Inches  Thickness 
Group;  Extrapolated  Test  Results. 
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Time,  hrs. 

igure  14(a).  Crack  Length  aft)  Versus  Time  t  for  11/32  Inches 
Thickness  Group;  Statistical  Model. 


Time,  hrs. 

Figure  14(b).  Crack  Length  aft)  Versus  Time  t  for  11/32  Inches 
Thickness  Group;  Extrapolated  Test  Results. 
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Figure  15(a).  Crack  Length  a(t)  Versus  Time  t  for  15/32  Inches 
Thickness  Group;  Extrapolated  Test  Results. 
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Figure  15(b).  Crack  Length  a ( t )  Versus  Time  t  for  15/32  Inches  Thickness 
Group;  Extrapolated  Test'  Results. 
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Figure  16(a).  Crack  Length  a(t)  Versus  Time  t  for  19/32  Inches 
Thickness  Group;  Statistical  Model. 


Figure  16(b).  Crack  Length  a(t)  Versus  Time  t  for  19/32  Inches 
Thickness  Group;  Extrapolated  Test  Results. 
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:igure  17(a).  Crack  Length  a(t)  Versus  Time  t  for  23/32  Inches 
Thickness  Group;  Statistical  Model. 


Time,  hrs 

Figure  17(b).  Crack  Length  a(t)  Versus  Time  t  for  23/32  Inches 
Thickness  Group;  Extrapolated  Test  Results. 
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Figure  18.  Correlation  for  Distribution  of  Time  to  Reach  Given  Crack 
Length  for  7/32  Inches  Thickness  Group;  Solid  Curves  for 
Statistical  Model  and  o  □  •  A  for  Extrapolated  Test 
Results. 
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TIME,  HOUR 

Figure  19.  Correlation  for  Distribution  of  Time  to  Reach  Given  Crack 
Length  for  11/32  Inches  Thickness  Group;  Solid  Curves  for 
Statistical  Model  and  o  □  •  a  for  Extrapolated  Test 
Results. 
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DISTRIBUTION  FUNCTION 
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TIME,  HOUR 


Figure  20.  Correlation  for  Distribution  of  Time  to  Reach  Given  Crack 
Length  for  15/32  Inches  Thickness  Group;  Solid  Curves  for 
Statistical  Model  and  o  □  •  &  for  Extrapolated  Test 
Results. 
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DISTRIBUTION  FUNCTION 
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Figure  21.  Correlation  for  Distribution  of  Time  to  Reach  Given  Crack 
Length  for  19/32  Inches  Thickness  Group;  Solid  Curves  for 
Statistical  Model  and  o  □  •  a  for  Extrapolated  Test 
Results. 
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Figure  22.  Correlation  for  Distribution  of  Time  to  Reach  Given  Crack 
Length  for  23/32  Inches  Thickness  Group;  Solid  Curves  for 
Statistical  Model  and  o  □  •  A  for  Extrapolated  Test 
Results. 
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The  crack  exceedance  curves  at  different  time  instants,  t,  are 
obtained  from  Figures  T 3( a) -17( a)  and  plotted  in  Figures  23-27  as  solid 
curves  for  the  statistical  model.  For  instance,  a  square  shown  on  the 
solid  curve  in  Figure  23  indicates  that  at  t=l . 5  hours,  the  probability 
that  the  crack  length  will  exceed  0.8  inches  (20.3  mm)  (abscissa)  is  0.4 
(ordinate).  The  corresponding  extrapolated  test  results  are  obtained 
from  Figures  13(b)-17(b)  by  drawing  a  vertical  line  through  the  time  of 
interest  and  the  results  are  displayed  in  Figures  23-27  as  circles.  As 
expected,  the  exceedance  probability  reduces  as  the  crack  length 
increases.  Again,  the  correlation  for  the  crack  exceedance  curves 
between  the  statistical  model  and  the  extrapolated  test  results  is  very 
reasonable. 
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Figure  23.  Probability  of  Crack  Exceedance  at  t=l . 5  Hours  for  7/32 
Inches  Thickness  Group;  Solid  Curve  for  Statistical  Model 
and  0  for  x  Extrapolated  Test  Results. 
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CRACK  LENGTH,  INCH 


Figure  25.  Probability  of  Crack  Exceedance  at  t=1.5  Hours  for  15/32 
Inches  Thickness  Group;  Solid  Curve  for  Statistical  Model 
and  0  for  Extrapolated  Test  Results. 
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Figure  26.  Probability  of  Crack  Exceedance  at  t=3.0  Hours  for  19/32 
Inches  Thickness  Group;  Solid  Curve  for  Statistical  Model 
and  0  for  Extrapolated  Test  Results. 
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CRACK  LENGTH,  INCH 

Figure  27.  Probability  of  Crack  Exceedance  at  t=2.5  Hours  for  23/32 
Inches  Thickness  Groujb ;  Solid  Curve  for  Statistical  Model 
and  0  for  Extrapolated  Test  Results. 
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SECTION  IV 
CONCLUSIONS 

A  fracture  mechanics-based  statistical  model  for  the  crack  growth 
behavior  of  engine  materials  under  sustained  loads  at  elevated  temperature 
has  been  proposed  and  applied  to  the  test  results  of  IN100,  a  superalloy 
used  in  F100  engine  disks.  The  method  of  maximum  likelihood  has  been 
employed  to  calibrate  the  fracture  mechanics  parameters  and  the  statistics 
of  the  model  from  the  test  results  of  the  crack  growth  rate  alone  without 
using  the  data  of  crack  length  versus  time.  A  correlation  study  is  per¬ 
formed  for  the  distribution  of  the  crack  length  as  a  function  of  time, 
the  distribution  of  time  to  reach  any  crack  length,  and  the  distribution 
of  the  crack  length  at  any  time.  A  good  correlation  has  been  demonstrated 
between  the  statistical  model  and  the  test  results. 

In  view  of  the  limited  experimental  test  results  presently  available, 
the  present  statistical  model  is  quite  practical.  This  is  because  the 
maximum  likelihood  estimates  of  a  few  fracture  mechanics  parameters  and 
the  statistics  of  the  model  require  neither  a  large  number  of  test  speci¬ 
mens  nor  a  homogeneous  data  base  for  crack  length  versus  time. 
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